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ABSTRACT 

Accretion among macroscopic bodies of ^km size or larger is enhanced significantly due to gravita- 
tional focusing. Two regimes can be distinguished. Initially, the system experiences runaway growth, 
in which the gravitational focusing factors increase, and bodies at the high-mass tail of the distri- 
bution grow fastest. However, at some point the runaway body dynamically heats its environment, 
gravitational focusing factors decrease, and runaway growth passes into oligarchic growth. Based on 
the results of recent simulations, we reconsider the runaway growth-oligarchy transition. In contrast 
to oligarchy, we find that runaway growth cann ot be approximated with a two component model (of 
small and large bodies) and that the criterion of llda fc Makinol ([1993), which is frequently adopted as 
the start of oligarchy, is not a sufficient condition to signify the transition. Instead, we propose a new 
criterion based on timescale arguments. We then find a larger value for the runaway growth-oligarchy 
transition: from several hundreds of km in the inner disk regions up to ~10^ km for the outer disk. 
These findings are consistent with the view that runaway growth has been responsible for the size 
distribution of the present day Kuiper belt objects. Our finding furthermore outline the proper initial 
conditions at the start of the oligarchy stage. 

Subject headings: planets and satellites: formation — Kuiper belt: general — protoplanetary disks — 
methods: statistical 



1. INTRODUCTION 

In planet formation theory, the process through which 
~km or larger bodies are formed is still not well under- 
stood, s ince a straightforward formation mechanism i s 
lacking (|Blum fc WurmI [20081 : IChiang fc YCTodhH [20091) . 
However, once these bodies appear gravity takes over 
the accretion process, fulfilling both the role as a stick- 
ing agent as well as an important accelerator for growth: 
the gravitational cross section between two bodies can 
become much larger than their geometrical cross section. 
This phenomenon is better known as gravitational fo- 
cusing and the gravitational enhancement factor in this 
regime is given by ^(wcsc/w)^, where Vcsc is the (mu- 
tual) escape velocity and v the velocity dispersion of 
bodies. Initially, gravitational focusing factors for the 
biggest bodies increase rapidly, and the system experi- 
ences runaway growth: Ucsc increases much faster than v. 
However, at a certain point the stirring capabilities of the 
biggest body will cause v to increase rapidly. Gravita- 
tional focusing factors decrease and the growth becomes 
self-regulated since the stirring rate of small bodies is de- 
termined by the same big body that accre tes them. Run- 
away growth has passed into oligarchy ^Kokubo fc Idal 

\T9m . 



Ilda fc Makinol (|1993[ ) have argued that the runaway 
growth-oligarchy transition takes place at the point 
where the stirring power of big bodies first exceeds that 
of the small bodies, i.e.. 



2EmM > 



(1) 



where Em is the surface density of big bodies of mass 
M and E^ that of the small bodies. Equation ([T]) 
can be transformed into a radius, i?ig/oiij indicating 
the turnover from runaway growth into oligarchy (see 
below, eq. [4]). Many works have adopted equation 
m as the start of t heir oligarchic calculations (e.g. , 
Thommes et"an[200l Ilda fc Linl[200l IChambersl[2006 . 



20081: iFortier et 3.1.1120071: iBrunini fc Benveriut'or boOS 



Miguel fc Bruninill2008l: iMordasini et al.ll2009D . 
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Ilda fc Makinol ([19931 ) and present a new expression 
for -Rrg/oii (eq. [13|). In runaway growth the col- 
umn density spectrum evolves into a power law, 

N(m) = (l/m)dJl/ dm o c m~^, where jo ~ —2.5 

(IWetherill fc Stcwafl [l99l IKokubo fc Idal [l99l [2OOOI : 
iBarnes et al.ll2009l . Fig.[TJi). However, for a, p > —3 in- 
dex, the stirring power lies at the high-mass end of the 
population; that is, during runaway growth the stirring 
power is already moving away from the initial mass (jtiq). 
Despite the stirring, the system continues in its runaway 
(fast) growth mode. The point is that equation ([Ij im- 
plicitly assumes that stirring occurs fast and outpaces 
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the accretion, which is true for ohgarchy but not for run- 
away growth. In other words, equation ([1]) is a necessary 
condition for ohgarchy but not a sufficient one. Instead, 
we will argue that the condition for the start of oligarchy 
is met when the stirring timescale in the two component 
approximation (Tvs, see below) drops below the accretion 
timescale that characterizes the runaway regime, T^g. 

In §[2] we first introduce key definitions and obtain the 
stirring timescale Tvs and the accretion timescale Tac for 
a two component system. In § [3] we present the results of 
our runaway growth simulations in terms of Ti-g, the ac- 
cretion timescale during runaway growth, which follows 
from our simulations. § 0] then presents the new tran- 
sition radius i?tr by equating the timescale expressions. 
We discusses a few implications and summarize in § [5] 

2. KEY DEFINITIONS AND TIMESCALES 

The Hill radius Rh and Hill velocity Vh of a single body 
of mass M and radius R are given by 



Rh = a 



f M 



V 3M* 



1/3 



Vh = Rh^, 



(2) 



with a the disk radius (semi major axis), M^, the mass 
of the star, and Q the orbital frequency at semi major 
axis a. The escape velocity is defined Wosc = \/2GM/R 
where G is Newton's gravitational constant. Using these 
definitions one can show that v^^^ = 6vl/a, where the 
dimensionless a is given by (cf. IGoldreich et aI1[200l 



R 

a= — = 7.5 X 10 
Rh 



Vau/ 



p 



g cm 



Mo 



(3) 

where we used M = AnpR^ /3 with p the internal density 
of the bodies. When discussing interactions one should 
use combined masses and radii in the above definitions 
(i.e., M = Mi+ M2 and R = Ri + R2) but this leaves 
equation ^ and the relation v^^^ = Gvf^/a unaffected. 

In the following we assume a two component model 
where a single big body of mass M and radius R interacts 
with smaller bodies of radius Rq and mass rriQ that dom- 
inate the solids surface density, E « Timo ^ Sm- The 
velocity dispersion v of the smaller bodies also dominates 
the relative velocity and we assume that Vh < v < Vcsc, 
where Vh and Vcsc correspond to, respectively, the Hill 
and escape velocity of the big body. This is the dispersion 
dominated regime and these assumptions are typical for 
oligarchy. Using Em = M/(27raAast) with Aost = ARh 
the width of the heating region and A an order-of-unity 
factor, one finds that the runaway growth-oligarchy tran- 
sition according to e quation ([T]) lies at a radius of (cf. 



sition according to ec 
iThommes erall[200l 



-^rg/oli " 



-,1/5 



4pa 
E 



= 94 km ( ^ 



1/5 



P 



10 g cm 



1/5 



/ a \2/5 



VAU J 



Ro 
10 km 



g cm 

3/5 



-2/15 



(4) 



In the dispersion-dominated regime i nclinations i ar e 
related to eccentricities e as i w e/2 (jlda et al.l I1993D . 
The accretion rate of the big bodies then becomes 



dM 
dt 



CttE- 



, 2GMR 



(5) 



(jlda fc Makinol [1993D where e is the rms-eccentricity of 
the small bodies, which is related to u as t; = eaO. Fur- 
thermore, C is a factor of order unity that takes into 
account the increased accretion rate under a distribution 
of velocities. With the above definitions we transform 
equation ([5]) into Hill units 
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dt a 

and define the accretion timescale in the 2-component 
approximation as 



Tac — 



1 dM 

M~dr 
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where we defined the dimensionless iCac = 2/9C. 

The time scale for viscous stir ring of the small bodies 
is given by (|Ida fc Makind[l993h : 



AttG^umM^ In A' 



(8) 



where In A is a Coulomb factor and hm the number den- 
sity of perturbers. This latter quantity is obtained from 
the assumption that each small bodies 'sees' one oligarch. 
The volume traversed by the small bodies is the product 
of 27ra (their circumference) , 2v/fl (the width in the ra- 
dial direction), and 2vz/^ (the vertical excursion). Here, 
Vz w v/2 represents the corresponding vertical veloci- 
ties as given by the inclinations i (=e/2) of the bod- 
ies. Thus, the effective number density of the single oli- 
garch is um = l/{2TTa){2v/il){2vz/^). Then, equation 
(l8l) transforms into 



with K,, 
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3. THE TIMESCALE FOR RUNAWAY GROWTH, Trq 

Figure [1] presents an example of the runaway 
growth parameter study we have recently conducted 
(jOrmel et al.l IMot ). These are the results of statis- 
tical simulations that include key physical processes 
like dynamical friction, viscous stirring, gas drag, 
and resolve t he serni-maj o r axis of the b odies (cf. 
Wcidenschillin g et al.l 119971 : iBromlev fc Kenv on 2006). 
In Fig. HJi the mass spectrum is shown at several times 
during the simulation run. After the radius of the most 
massive body, Ri{t), has increased by a factor 2 a new 
curve is plotted; the corresponding times are indicated 
in the legend. It is seen that the number density spec- 
trum N{m) evolves into a power law, N{m) cx to^p with 
p K, —2.5. Near the end of the simulation bodies have 
separated: the oligarchs. 

In Fig. [T)d additional information for this run is pre- 
sented as function of the evolutionary parameter Ri{t). 
The black curve shows the ratio of the most massive par- 
ticle (Ml) and m*, which is the mass- weighed average 
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Figure 1. (left) The mass spectrum dT,/dR for a simulation with gas drag conducted at 6 AU. After every factor of two increase of the 
radius of the most massive particle (-Ri) the distribution is plotted. The corresponding times are indicated in the legend. The dotted 
line in the left corner indicates a p = —2.5 power-law slope for the column density spectrum, N{m). The dotted line in the bottom-right 
corresponds to a distribution that contains one body per bin. The simulation is continued until Ri = 2 X 10'^ km. (right) Indicators of 
runaway growth: the mass of the most massive particle over the characteristic mass {solid black curve) and the maximum velocity in the 
distribution (vx) over the Hill velocity of the biggest body (vf^), which is a measure of the gravitational focusing {dashed black curve). 
These are plotted as function of the evolutionary parameter, Ri{t). Time is indicated by the grey curve on the second y-axis. The runaway 
growth timescale T^g is obtained from the linear portion of the curve, before it steepens at Ri = Rt^. 



of the distribution, = J m'^N{m)dm/ J mN(m)dm. 
For a two component model Mjjm^ shoul d increase if the 
syste m is in runaway growth (fLe e 2000; O rmel fc SpaansI 
|2008| ). It can be seen, however, that this ratio flattens 
near i?i ~ 200 km and decreases after ^600 km. The 
fluctuations in this curve are caused by merging of bod- 
ies of similar size. Also plotted is Vx/vh, the ratio of the 
largest random velocity of the smallest bodies over the 
Hill velocity of i?i . It first decreases (gravitational focus- 
ing increases) until Ri ~ 600 km, where Vx/vh reaches a 
minimum. This radius is denoted the transition radius, 
i?tr, indicated by the vertical dotted line. Time is plot- 
ted by the grey curve on the second, linear, y-axis. Af- 
ter Vx/vfi has reached its minimum this curve noticeably 
steepens. Only from this point onwards does dynamical 
heating win over the increase of Vh through accretion, 
and the accretion timescale increases. Bodies in other 
zones, spatially separated from Ri, then find the chance 
to catch up. 

The simulations show that during the runaway growth 
stage Ml (t) grows exponentially with time (i?i (t) is lin- 
ear on a semilog plot) and an exponential fit is appro- 
priate, Mi{t) oc exp(t/Trg) with T^g the runaway growth 
timescale (see Fig.[Tj), thin solid line). In Fig. [2] this 
quantity has been plotted for several runs performed at 
1, 6, and 35 AU. We express T^g in terms of the fiducial 
timescale irun — pRofS^ and derive the dimensionless 

Krg, 
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During the initial evolution T^g ^ Tac (eq. [7]) and also 
Trg ^ Tvs (eq. [9]). We find that the size and velocity 
spectrum that develops during runaway growth enables 
the biggest body to accrete particl es from all masses , 
i.e., also intermediate mass bodies (jOrmel et al.|[20Tol ). 
which decreases T^g compared to its 2-component esti- 
mate. Similarly, despite the fact that random velocities 
(vx) increase during runaway growth, the stirring by a 



Figure 2. (x-axis) The timescale for accretion during runaway 
growth, Trg for runaway growth simulations conducted at 1, 6, 
and 35 AU, and for models with and without gas drag and with 
(a simple prescription for) fragmentation, (y-axis) The dimen- 
sionless -R'rg, which gives Trg normalized to the fiducial time 
trun = ^op/Sfi, is much less than unity. Runaway growth is fast. 

single biggest body is simply insignificant to shape the 
dynamical evolution of the smaller bodies. Instead it is 
the ensemble of large and intermediate mass bodies that 
produces the stirring. A two component approximation 
is just too simple a picture for the runaway growth phase. 

4. THE RUNAWAY GROWTH/OLIGARCHY TRANSITION 

Although initially Tvs S> T^g and Tac ^ ^rg, during the 
runaway growth phase Vx/vh decreases and both Tac and 
(especially) Tvs rapidly converge on Tj-g. At the point 
where Tvs ^ Tj-g the runaway body starts to dynamically 
heat its environment at a rate faster than its previous ac- 
cretion rate. Thus, from this point the runaway growth 
can no longer 'outpace' the stirring. A single body dom- 
inates its neighborhood in terms of both velocity evolu- 
tion and accretion, meaning a much slower growth rate 
because accretion now becomes self-regulated. This is 
of course the key characteristic of oligarchy. Conversely, 
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Figure 3. The accretion timescale for the two component as- 
sumption Tac and Tvs (Eqs. [7] and [9]) evaluated at the point 
where Vx/vh reaches its minimum value, i.e., at Ri = Rtr and 
^x/vh = {v/vf^)ti, normalized to Trg. 

we may suspect that the condition Tvs — Tg heralds the 
end of the runaway growth stage and the transition to 
ohgarchy. 

In Fig. |3]we have plotted the ratio of these timescales 
over Trg at the point where Vx/vh is at its minimum (i.e., 
at Ri(t) = Rtr)- In evaluating T^^ and Tvs we have 
taken C = 3 (jGreenzweig & Lissauer||1992t K^c = 0.07) 
and InA = 3 (Stewart & Ida 2000, i^vs = 0.04). Because 
the Ri vs. v/vh curve (Fig. [T]d) differs from simulation to 
simulation (e.g., when including gas drag and fragmen- 
tation) the timescales Tac and Tev also evolve differently. 
However, from Fig. [3] it is clear that at Ri{t) — Rn- all 
the relevant timescales are similar: Tvs ^ Trg ~ Tac irre- 
spective of the simulation. Thus, the 2-component model 
becomes first applicable at the point where gravitational 
focusing reaches its maximum. This is the start of oli- 
garchy. From now on, Tj-g loses its meaning (runaway 
growth ceases), v^/vh increases, and accretion timescales 
are given by Tac (which also increases). 

Equating the relevant timescales, i.e., equations ([7]), 
and ([TU]), we solve for both i?tr and {vx/vii)tr 
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Equation ((TT|) supersedes -Rrg/oii of equation ([4]) as the 
new criterion between the runaway growth and oligarchy 
accretion phases. Using the above expressions for Kg^^ 
and Kys and inserting equation ([3]) for a we find that 
equation ()11|) transforms into 
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,(13) 



where the prefactor is ~3 times larger than equation 
(|H) . Any order-of- unity uncertainty in the relation T^g = 
Tac = Tvs can be regarded as an uncertainty in the K- 
factors, but this affects the prefactor only marginally. 
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Figure 4. The runaway / oligarchy transition, Rtv, against disk ra- 
dius according to equation l(4)l (Ida & Makino 1993, dashed lines) 
and equation I I13I1 (this work, solid lines). Plotted are _Rtr for 
a minimum mass solar nebula (MMSN) profile of surface den- 
sities {black lines) and a 10 times enhancement of S over the 
MMSN value {g rey lines). The MMSN profile is adopted from 
INakagawa et all l| 198(|) . 

More fundamentally, equation ([T^ shows different de- 
pendencies on i?o, S, and, especially, a. If i?o is large, 
e.g., i?o ^ 500 km, i?rg/oii and i?tr give approximately the 
same values at 1 AU. A large Rq is the favored outcome of 
recent numerical simulation s involving a turbulent laye r 
of densely packed boulders (jJohansen et al.ir2007ll2009l) . 
Then, the domain of runaway growth is rather limited 
(a factor three in size); but it still produces ~10'^ km 
objects. 

In Fig. [4] we illustrat e the run a .way g rowth/oligarchy 
turnover point of Ida fc Makinol (| 1993( 1 (eq. [Ij , dashed 
lines) and equation l|13p (solid lines) as function of dis- 
tance a for an initial planetesimal radius of Rq = 10 km 
and an internal density of p = 3 g cm~^. For E(a) we 
adopt a Minimum Mass Solar Nebula (MMSN) profile 
(^cidcnschilling 1977) (black lines) and a profile ten 
times larger than this (grey lines). In Figure |4] it can 
be seen that the difference between the two criteria be- 
comes largest in the outer disk (a ^ 10 AU): whereas 
according to the il da fc Makinol (|1993[ ) prescription the 
transition takes place at radii of ~100 km, the new cri- 
terion shifts it upwards to '^10'^ km. 

5. DISCUSSION 

According to our findings, the transition between the 
runaway growth and oligarchic growth phases is charac- 
terized by the following properties: 

• A power-law size distribution of mass index p « 
—2.5, extending from the initial size Rq to the tran- 
sition size Rtr as given by equation ()13|) : 

• The random velocity v of the planetesimal bodies 
(of size ^Rq), via equation (|12|): 

• The timescale, via equation (fTO|) . Here, Tj-g must 
be multiplied by a term ~log(i?tr/^o) to account 
for the several e-foldings of enjoyed growth. Then 
we obtain t = Ttr ~ pRo/T,o^l for the time un- 
til the transition. It is remarkable that this short 
timescale depends on the initial conditions only, a 
result that is unique to the runaway growth phase. 
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We assess the imphcations of these findings for the 
broader context of planet formation. First, we do 
not expect the final timescales of core formation to 
be much influenced by the new transition radius, since 
these are set by the much slower oligarchy stage that 
supersedes runaway growth. Indeed, semi-analytical 
studies of oligarchic growth required additional mech- 
anisms like gas damping, fragmentation, or migra- 
tion to produce embryos on reasonable times cales (e.g., 
iBrunini fc Benvenutol 120081 : 'Chambers' '20081). On the 
other hand, gap formation (Rafikov 2003) will increase 
formation timescales^ 

Recently, iLevison et al.l (|201C1[ ) investigate core- 
formation scenarios using A'^-body techniques. The aim 
of that study was to understand how efficient planetesi- 
mal accretion proceeds during the phase where an Earth- 
size planetary embryo has to grow to a mass of ~10Mgg 
to be able to accrete the nebula gas and become a 
gas giant. They recognize the importance of processes 
like planetesimal scattering, planetesimal orbital decay 
due to gas drag (both processes quench the growth), 
and planetesimal-driven embryo migration (which is con- 
ducive to growth). The outcome of these simulations, 
furthermore, is found to depend on the initial setup of the 
simulation, i.e., the size distribution of the planetesimals. 
Given the significance that is attached to the wholesale 
redistribution of matter, it would also be of interest to 
assess the importance of these effects for the early oli- 
garchy stage, e.g., to perform A^— body simulations with 
embryos radii starting at the transition mass Rt^- 

The second implication of our study concerns the 
Kuiper belt. Our results strongly suggest that the Kuiper 
Belt is primarily the product of the runaway growth 
phase. First, assuming that the initial surface density 
is approximately MMSN or larger (S ~ 0.1 g cm^^), the 
biggest ^10^ km bodies (plutinos) can be produced by 
runaway growth. Second, the observed mass distribu- 
tion N{m) for the largest Kuiper belt objects (KBOs) 
obeys a power-law with p w —2.5. Recent studies find a 
power-law size index q (as in N{R) oc i?~*) of 4.8 ± 0.3 
{p = -2.3 ± 0.1; iFraser fc Kavelaarsll2009l ) and 4.5l:J ° 

(p = -2.2toi; iFuentes fc Holma n' '2008^ . (The size dis- 
tribution of the Kuiper belt's smallest bodies is collision- 
ally dominated, though, and q is much lower.) The large 
q, together with the fact that the KBO size distribution 
is continuous rather than bimodal as in the end stage of 
Fig. argue that it has evolved only very little since 
runaway growth and has not been significantly shaped 
by oligarchic growth. 

In order to produce the KBOs in a sufficiently 
short time span, Kuiper Belt formation s cenarios (e.g., 
IKenvon fc Luul 119981: IChiang et all I2007D assume that 
the initial belt contained much more mass than the 
0.01 Mm (E ^ 0.0 01 g cm-2) that is present today 
(jBernstein et al.ll2004[ ). Neptune formation and/or an- 
other dynamical shakeup event (as in t he Nice model) 
subsequently depleted 99 % of its mass (iFord fc Chiand 
[20071 ILevison et al.l[200l . Using our results, we can ver- 
ify these findings. Assuming a = 35 AU, we find 



which readily shows the need for an enhanced surface 
density over the current one: for E ~ 0.001 g cm~^ ei- 
ther Rtr becomes too low or Ttr too long. Thus, we 
conclude that the KBO size distribution as seen today 
is consistent with a scenario of being a leftover product 
of the initial runaway growth phase and has since been 
de pleted. These findin g s are i n line with the simulations 
of Kenvon fc BromlevI (|2008l I2OO90 . where growth also 
stalls after ^10^ km bodies have been formed. Further 
growth is impeded since oligarchic accretion timescales 
become too long, even at enhanced surface densities. 

We acknowledge the helpful comments of the (anony- 
mous) referee and his/her encouragement to put our re- 
sults into a broader context. C.W.O. is supported by a 
grant from the Alexander von Humboldt foundation. 
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